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Physics beyond the Standard Model predicts the possible existence of new particles that can 
be searched at the low energy frontier in the sub-eV range. The OSQAR photon regeneration 
experiment looks for “Light Shining through a Wall” from the quantum oscillation of optical photons 
into “Weakly Interacting Sub-eV Particles”, such as axion or Axion-Like Particles (ALPs), in a 9T 
transverse magnetic field over the unprecedented length of 2 x 14.3 m. In 2014, this experiment has 
been run with an outstanding sensitivity, using an 18.5 W continuous wave laser emitting in the 
green at the single wavelength of 532 nm. No regenerated photons have been detected after the wall, 
pushing the limits for the existence of axions and ALPs down to an unprecedented level for such a 
type of laboratory experiment. The di-photon couplings of possible pseudoscalar and scalar ALPs 
can be constrained in the nearly massless limit to be less than 3.5-10“® GeV“^ and 3.2-10“® GeV“^, 
respectively, at 95% Confidence Level. 
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I. INTRODUCTION 

Possible extensions of the Standard Model (SM) of par¬ 
ticle physics are not restricted to the high energy frontier. 
There is also a growing interest for the search of Weakly 
Interacting Sub-eV Particles (WISPs), with much weaker 
interactions and masses below the eV. One emblematic 
example is the axion, a pseudoscalar boson arising from 
the spontaneous breaking of a global chiral symmetry 
1/(1)A, postulated to dynamically solve the strong GP 
problem [1-3]. Axions and Axion-Like Particles (ALPs) 
are predicted in supersymmetric theories [4], in string 
theory [5, 6], and in the Gonformal Standard Model [7]. 
ALPs can be scalar as well as pseudoscalar and theorized 
to couple to the SM through a variety of mechanisms, 
giving rise in particular to a two-photons vertex. WISPs 
also include light bosons of gauge groups under which the 
SM particles are not charged (hidden sectors), which may 
interact with the SM through gravity, kinetic mixing or 
higher order quantum processes [8] . The interest aroused 
by WISPs goes beyond particle physics. As earlier hy¬ 
pothesized for the axion [9] , they provide alternative can¬ 
didates for Dark Matter (DM) [10-12]. Moreover, they 
might explain a number of astrophysical puzzles, such as 
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the universe transparency to very high energy photons 
(> 100 GeV) [13], the anomalous white dwarf cooling [14] 
or the recently discovered gamma ray excesses in galaxy 
clusters [15]. In contrast to Weakly Interacting Massive 
Particles (WIMPs), which can be searched for at TeV 
colliders such as the Large Hadron Gollider (LHG) at 
GERN, the detection of WISPs requires to have recourse 
to dedicated low-energy experiments. Several method¬ 
ologies exploiting the existence of a di-photon coupling 
have been proposed based on lasers, microwave cavities, 
strong electromagnetic fields or torsion balances [16, 17]. 

The OSQAR (Optical Search for QED Vacuum Bire¬ 
fringence, Axions and Photon Regeneration) experiment 
at GERN, is at the forefront of this low-energy frontier 
of particle/astroparticle physics. It combines the simul¬ 
taneous use of high magnetic fields with laser beams in 
distinct experiments. One of its setups uses the “Light 
Shining Through a Wall” (LSW) method for the search 
of the ALPs [18]. A pioneering work in this line excluded 
ALPs with a di-photon coupling constant larger 

than 6.7 • 10“^GeV“^ for masses below 10“^ eV [19]. 
These exclusion limits were later extended by other LSW 
experiments to > b.5 • 10“^ GeV“^ for masses be¬ 
low 5 X 10“^ eV [20] and more recently tightened to 
OA'y'y > 5.7-10“^ GeV“^ for masses below 2x10“^ eV [21]. 
New results are here reported, obtained from the 2014 
data-taking campaign of the OSQAR LSW experiment. 
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Compared to the previous experimental setups of OS- 
QAR [ 21 , 22 ], a more powerful laser source and a detector 
with higher sensitivity have been used. In addition, an 
improved data analysis strategy has been implemented. 

II. EXPERIMENTAL SETUP AND DATA 
TAKING 

LSW experiments are based on the combination of 
two factors. The first is the transparency of the WISPs 
to photons barriers, owing to the weakness of the in¬ 
teractions of the WISPs with the particles of the SM. 
The second is the photon-to-WISP and WISP-to-photon 
quantum oscillation, which would arise from the interac¬ 
tions of WISPs with photons [18, 23, 24]. In the case of 
the sub-eV ALPs, the method benefits from their two- 
photon vertex, inducing oscillations with optical photons 
in a transverse magnetic field. Such a field can be rep¬ 
resented as a sea of virtual photons, whose interaction 
with real photons (resp. ALPs) can produce real ALPs 
(resp. photons). The mechanism is similar to the Pri¬ 
makov process of the production of neutral mesons by 
high-energy photons in a strong electric field [25]. The 
effective Lagrangian density of the interaction of an axion 
or a pseudoscalar ALP (PS-ALP) field A with the elec¬ 
tromagnetic field is written generically in the form 

i^int =-^QA-f-f A A E ■ B (1) 

where F^*'' = is the dual of and the 

constant ^a 77 stands for the effective axion or ALP di¬ 
photon coupling. With a Scalar ALP (S-ALP) field A 
the interaction with the field takes the generic form 

Ant = -^^A77 a F^j.F^'^ = ^A77 A -{E^ - B^) (2) 

Accordingly, either PS-ALPs or S-ALPs could potentially 
be created when a beam of linearly polarized photons 
propagates in a transverse magnetic field 5, depend¬ 
ing on whether the polarization is parallel to the mag¬ 
netic field {Ej II B) or perpendicular {E^ ± B). If 
an optical barrier is placed downstream to the beam, 
all unconverted photons will be absorbed while ALPs 
would traverse the optical barrier. By applying a sec¬ 
ond magnetic field in the regeneration domain beyond 
the wall, the inverse Primakov process can convert the 
ALPs back into photons, which can be subsequently de¬ 
tected (Fig. 1 ). The probability of an ALP-to-photon 
{A 7 ) or of a photon-to-ALP (7 ^ A) conversion is 
given by [18, 23, 24]: 

P-f^A = sin (3) 

in units of Heaviside-Lorentz system {h = c = 1). 
q = \k^ — kA\ stands for the momentum transfer, where 


Optical barrier 



FIG. 1. Principle of the OSQAR LSW experiment. 


k^ = uj is the momentum of the photon of energy uj and 
kA = a/(cj^ — m\) the momentum of the ALP of mass 
rriA- The overall probability of the photon regeneration 
is Pj^A^'y = 1 as arising from the two consecu¬ 

tive conversions: 7 ^ A followed by A ^ 7 . The flux of 
the regenerated photons to detect is then given by 

dN ^ X 2 / N 

-^ = -1 (A^a) (4) 

where P represents the power of the incoming photon 
beam and 7 stands for the photon detection efficiency. 
dN/dt is proportional to the 4-th power of the field inte¬ 
gral BL^ whence the need for strongest magnetic field B 
over longest optical path L to achieve highest sensitivity. 

The OSQAR LSW experiment uses two spare LHC 
dipole magnets, the first for the ALPs production and 
the second for the photon regeneration. Each magnet 
is cooled down to 1.9 K with superfiuid He and provides 
an uniform transverse magnetic field with a strength of 
9 T over a magnetic length of 14.3 m, thus giving rise to 
a magnetic field integral of BL = 128.7 Tm. The aper¬ 
ture of both magnets have been pumped, using turbo- 
molecular pumping groups, down to 10 “^ — 10 “^ mbar. 

A continuous wave optical power of 18.5 W at a single 
wavelength of 532 nm (2.33 eV) was delivered by a diode- 
pumped solid-state laser (Verdi V18 from COHERENT 
Inc.). A beam expander telescope has been used to re¬ 
duce the laser beam divergence. The photon beam is lin¬ 
early polarized with a vertical orientation parallel to the 
magnetic field for the search of PS-ALPs. A A /2 wave- 
plate with antireflective (AR) coating layers was inserted 
at the laser exit in order to align the photon polarization 
in the perpendicular direction for the search of S-ALPs. 

The laser beam at the exit of the second magnet aper¬ 
ture was focused by an optical lens on a thermoelectric 
cooled Charge-Coupled Device (CCD) with an AR coated 
window (DU934P-BEX2-DD from ANDOR Co.). The 
CCD chip has an active area of 13.3 x 13.3 mm^ and 
is composed by a 2D array of 1024 x 1024 square pix¬ 
els of size 13 X 13 The laser beam was focussed 

at 95% on an area covering not more than 4 pixels. 
The CCD was cooled down to the temperature range 
[—95° C, —92°C]. The dark current is 0 . 0012 e“/pixel/s 
at —95° C and 0.0020 e“/pixel/s at —92° C. The read¬ 
out noise was determined at 2.5 e“rms/pixel from frames 
recorded in minimal time with closed shutters. These 
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measured CCD parameters are in agreement with those 
provided by the manufacturer. The overall photon 
detection efficiency was explicitly measured to 77 = 
0.56 ± 0.02 ADU/photon, taking into account the losses 
coming from optical elements including windows, plates, 
lens, CCD quantum efficiency of 0.88 at 532 nm and sen¬ 
sitivity of 1.3e“/Analog-to-Digital Unit (ADU). The to¬ 
tal optical power loss of the whole optical setup is thus 
approximately 17%, in good agreement with the charac¬ 
teristics of the various optical elements in use. 

The data taking for the ALPs search reported in this 
paper was performed in August 2014 and corresponds to 
a total of 119 experimental runs, as detailed in Table I. 


Search for 

E~i 

B 

P 

Nt 

Nr 

W(95%) 

PS-ALPs 

II B 

9T 

18.5W 

59 

18 

0.64 mHz 

S-ALPs 

T B 

9T 

18.5W 

60 

12 

0.45 mHz 


TABLE I. Summary of the experimental runs. stands 
for the electric component of the linearly polarized photon 
field, B the static magnetic field, P the laser power, Nt the 
number of recorded runs and Nr the number of rejected runs. 
W(95%) is the Bayesian threshold of non detection at 95%. 

Each experimental run was composed of two frames of 
5400 s exposure time separated by a 100 s pause. The 
recording time for a single frame was defined to optimize 
the signal over noise ratio, taking into account the filter¬ 
ing and removal procedure of cosmic ray signatures. Be¬ 
fore and after each run, the position of the laser beam was 
precisely measured on the CCD with the optical power 
strongly reduced. Starting with a laser beam power set¬ 
tled to 3 W, a variable beam splitter has been introduced 
to attenuate the beam typically by a factor of 1/500. 
To carefully check the stability of the laser beam spot, 
three frames were recorded each time with 0.01 s expo¬ 
sure separated by a 120 s pause. Subsequently the optical 
barrier was introduced and the attenuator was settled to 
its minimum before ramping up the laser power up to 
18.5 W. If the experimental conditions were not optimal 
or if cosmic rays have impacted the signal region then 
the experimental run was rejected (Table I). The signal 
region is defined as the pixels of the CCD where recon¬ 
verted photons are expected to be detected, taking into 
account slight displacements of the laser beam spot. 


III. DATA REDUCTION 

The detailed analysis of all beam positions during the 
full data taking period clearly confirmed that environ¬ 
mental temperature variations induce only minor low- 
frequency displacement/deformation of the CCD support 
and therefore only a small relative shift of the beam spot 
on the detector, typically less or equal to one pixel per 
hour. The beam spot displacement is unidirectional, with 
no oscillation observed after the exposure time of a full 
run, i.e. after 2 x 5400 s+100 s = 10900 s. The beam spot 
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EIG. 2. (Color online) Initial and final positions of the at¬ 
tenuated laser beam on the CCD for Run-90 devoted to the 
search of S-ALPs T B). The signal region covers 7 pixels. 


sizes before and after each experimental run are obtained 
from least squares fits with a two-dimensional Gaussian 
distribution. To define conservatively the signal region, 
twice the maximum value of the gaussian widths com¬ 
puted from the intitial (i) and final (/) reference frames 
is retained: 2crRun = 2 x maxfcr^, a/]. If the laser po¬ 
sitions before and after a run do not coincide, all pix¬ 
els within a 2crRun orthogonal distance along a straight 
line between the final and initial positions are added to 
the signal region. An illustrative example is provided in 
Fig. 2. It follows that the number S of pixels that defines 
the signal region can vary from run to run (Table II). A 
cumulated count Ni can be associated to each Run-i, by 
adding the counts collected during the exposure time of 
2 X 5400 s for ALPs detection on all the pixels of the sig¬ 
nal region, once this one is carefully determined. Ni is 
the relevant quantity that would reveal the regenerated 
photons if it is proved in excess over background noises. 
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TABLE 11. Distribution of sizes S in unit of pixel of the signal 
regions. Nrs (resp. Ns) counts the number of runs devoted 
to the search of PS-ALPs (resp. S-ALPs) for which the same 
given size S was experimentally deduced for the signal region. 

One of the advantages of a small and well-defined sig¬ 
nal region is the possibility to use the remaining pixels, 
which are not exposed to possible ALPs signals, to fully 
characterize the background without repeating dedicated 
background acquisition runs [21, 22]. These remaining 
pixels are clustered in size and shape that correspond to 
the signal region. Cumulated counts Uij are obtained for 
each cluster j in the background region of each Run-i, by 
adding the counts collected on all the pixels of the chosen 
cluster, similarly to the cumulated count Ni of the sig¬ 
nal region. The distribution of the Uij for a given Run-i 
then provides an estimate of the expected background 
contribution to Ni. As it is well known for CCDs, the 
width of this distribution arises mainly from two inde¬ 
pendent sources, namely dark-counts that increases with 
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exposure times and read-out noise that arises from the 
analog-digital conversion process of the readout-system. 

Before analyzing the cumulative counts, three correc¬ 
tions have been applied to each recorded frame pair of 
2 X 5400 s. The first one aims to correct long wavelength 
distortions of the spatial distribution of the background 
arising from the minor heterogenity in the cooling of the 
CCD chip. Corrections of these background drifts are 
achieved by median filtering. The objective of the sec¬ 
ond correction is to reduce the so-called fixed pattern 
noise evolving from possibly unequal pixel biases. As 
this noise manifests only during readout, and is inde¬ 
pendent of the exposure time, it has been estimated by 
capturing dark frames (z.e. shutter closed) with the min¬ 
imal exposure time (10“^ s in our case). Ten frames have 
been recorded under these conditions and have been aver¬ 
aged pixel-wise to define a bias frame, which is subtracted 
from every frame used in the data analyses. This correc¬ 
tion also subtracts the general constant offset of ~ 3300 
counts on each pixel imprinted by the chip in each frame 
of 5400 s. The last correction step removes those pixels 
which were hit by cosmic rays, using a short wavelength 
filtering process. Whenever a pixel count exceeds the 
significant threshold of the distribution of pixel counts in 
the background region this pixel is rejected together with 
its nearest neighbours. In case the signal region is hit by 
a cosmic ray, the overall frame is rejected. 

Once all correction steps have been applied to a given 
Run-i, the cumulated counts Ni in the signal region and 
Uij of the clusters j in the background region are evalu¬ 
ated. As illustrated in Fig. 3 for the Run-90, the distri¬ 
bution of the cumulated counts Uij for every Run-i is a 


0 



cumulated counts 


FIG. 3. (Color online) Distribution of the cumulated 
counts ngoj of background clusters j of size and shape identi¬ 
cal to the signal region shown Fig. 2 for Run-90, after ap¬ 
plying all background corrections (see text). Gaussian fit 
gives mean /XgQ® = 125.24 ± 0.05 and standard deviation 
cTgo® = 14.77±0.04 with reduced chi-squared x^/n.d./ = 1.04. 
The red line represents the cumulated count Ago = 119 of the 
signal region. It crosses the Gaussian distribution at about 
2000, which corresponds to the number of background clusters 
with the same cumulated count as the signal region. 


gaussian with mean and standard deviation in 
good agreement with the total noise of the CCD expected 
from dark-counts and read-out noise contributions. 


IV. BAYESIAN ANALYSIS 


A detection limit on the rate of reconverted photons 
dN/dt can be extracted from the gaussian width of cumu¬ 
lated background cluster counts over the pixel-wise sum 
of all the frames for a given photon polarization. This 
requires choosing a signal region of the same size for all 
the runs, which necessarily must be the largest one that 
is S' = 20 pixels according to Table II. A slightly bet¬ 
ter detection limit can be obtained by rejecting the few 
runs with too large size of signal region, but at the sta¬ 
tistical cost of the counts inherent to the corresponding 
decrease of the total exposure time. This strategy used 
in previously published analyses of similar experiments, 
e.g. Ref. [22], is not optimal. A detection limit at two 
times the standard deviation of the background fluctua¬ 
tions in the order of 2 mHz (2 x 10“^ADUs“^) can be 
obtained with this approach. This does not constitute 
a real improvement with respect to a single run with a 
smaller size of signal region such as Run-90 for which 
S = 7 pixels, because the increase in cluster size leads to 
larger background fluctuations for such a run. 

In order to optimize the detection limit by considering 
real signal regions of different sizes from run to run (Table 
II), a Bayesian approach has been implemented (cf. Ref. 
[26] pages 472-487). The following Likelihood model has 
been considered 


C oc 




with the rate of reconverted photons dN/dt as the sig¬ 
nal parameter. The index i iterates over the integer 
number of runs corresponding to a certain polariza¬ 
tion state and JV represents the Gaussian background 
parametrization of the i-th frame pair including an ad¬ 
ditional Poissonian signal contribution of 



FIG. 4. (Golor online) Posterior probability distribution 
function p(dN/dt \ Data) of reconverted photon rate dN/dt 
for incoming photon polarization parallel (E^ || B) (Left) and 
perpendicular {E^ T B) (Right) to the magnetic field B. 
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V. EXCLUSION LIMITS 

As a conclusion, no significant excess over the back¬ 
ground expectation has been observed in the signal re¬ 
gions neither for the parallel nor the perpendicular po¬ 
larization. Exclusion limits on the di-photon coupling 
strength gA^y^ and the ALPs mass tua have then been 
derived from Eqs. 3 & 4. 

The 95% Confidence Limit (C. L.) on the reconverted 
photon fiux is derived via the posterior distribution of 
the signal parameter dN/dt (Eig. 4). The results for the 
PS-ALPs and S-ALPs searches are summarized in Eig. 6, 
which also shows the exclusion limits reported by ALPS 
collaboration [20] . A limit for the di-photon couplings of 
^/A77 < 3.5 • 10-^GeV-^ and < 3.2 • 10“^ GeV-^ 
is obtained for the pseudoscalar and scalar searches for 
rriA < 2 X 10“^ eV, respectively. These results estab¬ 
lish the most stringent constraints on ALPs searches in 
the nearly massless limit obtained so far in LSW experi¬ 
ments. 


FIG. 5. (Color online) Posterior probability distribution 
function p{dN/dt \ Data) of reconverted photon rate dN/dt 
for incoming photon polarization parallel \\ B) (Top) and 
perpendicular T B) (Bottom) to the magnetic field 5, T 

with an artificially imposed fake-signal of 1 x 10“^ ADU s“^. ^ 

CD 

c^ 

expectation value dN/dt^ multiplied by the frame pair 
exposure time = 2 x 5400 s. The numerical inte¬ 
gration was performed via Markov-Ghain-Monte-Garlos 
using twenty million samples for each photon polariza¬ 
tion state, leading to a negligible uncertainty due to the 
integration process. The posterior probability distribu¬ 
tion function p{dN/dt \ Data) of reconverted photon rate 
dN/dt deduced from the considered likelihood model as¬ 
suming a flat prior distribution for the signal parameter 
dN/dt are displayed Eig. 4 for the two linear polarizations 
of the incoming photon beam. 

The consistency of the overall data analysis, including 
filtering and statistical methods was checked by impos¬ 
ing into all selected raw data frames fake-signals with a 
rate of 1 mHz corresponding to an hypothetical ALP with 
a di-photon coupling constant gA^y^ = 3.8 10“^ GeV“^ 
in the nearly massless limit {rriA ^0). A flux dN/dt = 
1.14±0.28mHz was inferred from the posterior probabil¬ 
ity distribution function p{dN/dt \ Data) of reconverted 
photon rate dN/dt for incoming photon beam linearly 
polarized parallel to the magnetic field and dN/dt = 

0.95 ± 0.33 mHz for incoming photon beam linearly po¬ 
larized perpendicular to the magnetic field (Eig. 5). 



FIG. 6. (Color online) Exclusion limits for the searches of PS- 
ALPs and S-ALPs at 95% C.L. obtained in vacuum by the 
present OSQAR LSW experiment. Also reported for compar¬ 
ison is the latest result of ALPS LSW experiment in vacuum. 
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